Long-term effect of inter-mode transitions in quantum Markovian process 
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We study the Markovian process of a multi-mode open system connecting with a non-equilibrium 
environment, which consists of several heat baths with different temperatures. As an illustration, 
we study the steady state of three linearly coupled harmonic oscillators in long time evolution, two 
of which contact with two independent bosonic heat baths with different temperatures respectively. 
We show that the inter-mode transitions mediated by the environment is responsible for the long 
time behaviour of the dynamics evolution, which is usually considered to take effect only in short 
time dynamics of the system immersed in a equilibrium heat bath with a single temperature. In 
non-equilibrium systems, the ignorance of the inter-mode transitions may give rise to some counter- 
intuitive result, thus they cannot be neglected. 
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I. INTRODUCTION 

When a small system contacts with a simple environ- 
ment, i.e., a heat bath in canonical equilibrium with a 
temperature T, it would approach its canonical state with 
the same temperature T. This dynamics process, which 
erases the information of the initial state, is called ther- 
malization [11-01 • The same thermalization process may 
also happen for a composite system, when all of its sub- 
systems are coupled with the same common heat bath. 
The composite system would approach its canonical state 
of temperature T with respect to its total eigen modes. 

However, nonequilibrium systems, e.g., the photosyn- 
thesis system, arc more general in nature, and exhibit 
more rich physics [6- 8] . A typical example is a composite 
system connecting with more than one heat baths with 
different temperatures. For a long-term evolution, the 
open composite system would not approach its canon- 
ical thermal state, but still it would be stabilized to a 
certain steady state. We call this process non-thermal 
stabilization. 

In an open system, both equilibrium and non- 
equilibrium case as we mentioned above, the energy ex- 
change with environment would mediate the transitions 
between the normal modes of the total system. It is usu- 
ally believed that this inter-mode transition only takes 
effect on the dynamics of transient evolution within the 
time scale determined by the time-energy uncertainty 
P, I9I412I] . If this transition oscillates too fast in the re- 
laxation time, averagely it has no contribution and can 
be neglected, which is known as rotating-wave approxi- 
mation (RWA). 

In this paper we consider whether or not such inter- 
mode transitions have a long-term effect on a composite 
open system contacting with a complicated environment, 
which consists of several independent heat baths with 
different temperatures. As an illustration, we study the 
steady state of three linearly coupled harmonic oscilla- 
tors (HOs), two of which contact with two independent 
bosonic heat baths with different temperatures respec- 



tively. Specially, we are interested in the steady state of 
each local HO. It is found that if the inter-mode transi- 
tion were ignored, there would be some counter-intuitive 
results in the long time steady state, especially when the 
frequency detuning and coupling strength of the local 
modes are small. In this regime, the coherent transitions 
among the total normal modes oscillate so slowly that 
they result in long-term effects, thus cannot be neglected. 

The paper is arranged as follows. In Sec. II, we setup 
the model of the coupled system and give a master equa- 
tion. In Sec. HI, we show the stabilization result and 
make some analytical discussion by eliminating the de- 
gree of freedom of the mediating data bus. We propose 
a possible implementation in Sec. IV. The calculation is 
assisted by some properties of the characteristic descrip- 
tion of Wigner function and Fokker-Planck equation. We 
leave these tricks in the appendices. Finally, summary is 
drawn in Sec. V. 



II. MODEL SETUP 

To study the long-term dynamics of a composite sys- 
tem coupled to a complicated environment, we use the 
coupled HOs system as an illustration. The system we 
study here is illustrated in Fig. 1. Two HO separate from 
each other for a distance. They are coupled indirectly by 
another mediating system, which is also modeled as a 
HO here. Effectively, we suppose the mediating HO does 
not contact with any environment. The three oscillators 
system can be described by a quadratic coupled Hamil- 
tonian Hs = Hq + V, where 



(1) 



and Hq describes the free Hamiltonian with local modes 
respectively defined by annihilation operators ol , or and 
b; V describes the coupling among the local modes. 

We assume the two oscillators locate remotely at differ- 
ent places, thus they may suffer from independent baths. 
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FIG. 1: (Color online) Demonstration of the coupled oscilla- 
tors system. Two remotely located HOs are indirectly coupled 
by another mediating one. The HOs at the two ends con- 
tact with independent heat baths with different temperatures 
Tl/r- 



We also assume each bath stays at a canonical thermal 
state with a temperature /r- The whole system can be 
described by the total Hamiltonian T-L = Hs + Hb + Vsb, 
where 



VsB^ air^ + a^Tl (2) 



a=L,R 



and To- = Ck„ . Hb is the free Hamiltonian of 

the two boson heat baths, each of which is modeled as a 
collection of boson modes, described by the boson annihi- 
lation operators Ck„. Vsb represents the linear coupling 
between the system and the environment. 

We need to derive a master equation to study the dy- 
namics of the open system. Actually for the coupled 
HO system, a correct treatment of the master equation 
should be based on the normal modes of Hs, but not the 
local modes a^/R and h. Otherwise, it may give rise to 
some counter-intuitive results. Thus, we diagonalize the 
Hamiltonian Hs as. 



= a^^ • n • a = 



gL 

gh i^m gR 
gn 




3 
i=l 



(3) 



where a = (a^, b, clr)'^ , and we also denote it as 
(cii, (22, 03)"^ hereafter (with redefined indices 1,2,3 for 
ol, &, CLR respectively). A — U ■ & — [Ax, A2, A^)'^ for 
Aj's being the normal modes. U-il-U^ = diag {ei, £2, £3} 
gives the eigen frequencies of the normal modes. Al- 
though the normal modes are decoupled from each other 
in the isolated Hs, we can see below that the environment 
could mediately induce some effective couplings between 
these normal modes. 

With the above notations, in Appendix A we derive a 
master equation to describe the long-term dynamics of 
the open system via Born-Markovian approximation ^ . 



In Schrodinger's picture, it reads as 
A 



where 



J2^i2^pA]^{A]A,pU) 



At, ~UaU;,[NL{ei) + NLie/)] 



(4) 



U^3U*[NR{e,) + NrIe,)], 



(5) 



=^UaU*[NL{ei) + NUe/) + 2] 



2 """^ 
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+ ^UaU*^[NR{e,) + NR{e,) + 2]. 



Here, 70- (^i) = 27rJo.(ei) characterizes the coupling 
strength with each bath, and Ja{^) = X^k l5k„| 5{ui — 
cjk^) is the coupling distribution. For the usual case, 
we can assume that 70- (w) ~ 7o- does not depend too 
much on lo and can be treated as constant. N„{uj) — 
[eyLY>{ijj / kTa) — is the Planck distribution for a ~ 
L,R. 

It is observed from the above master equation that the 
environment indeed induces an effective coupling between 
two normal modes Ai and Aj. Af, measure the transi- 
tions of the normal modes. This environment-mediating 
effect can be understood in the following way. The cou- 
pled system exchanges energy with environment through 
the interaction Vsb ■ Immediately after the normal mode 
Ai of the system emits an energy quanta Si to the envi- 
ronment, another process may happen in succession that 
the normal mode Aj absorbs back ej from the environ- 
ment. Also, there is another possibility that mode Ai 
first absorb energy from the environment and then the 
energy is emitted through mode Aj. Thus, with the me- 
diation of the environment, different normal modes Ai^s 
of Hs are coupled. In fact, this effect of environment me- 
diated coupling can be also found for two modes coupled 
to a common heat bath . For comparison, we can also 
understand our result according to the case of common 
bath. Although there are two independent baths, the 
two baths are combined into a single one by the coupling 
between the local modes. 

The transition terms with i j describes the effective 
coupling between different normal modes. In the interac- 
tion picture, these terms would contribute an oscillating 
factor exp[±j^eij t] resulted from the energy difference 
of the modes Ai and Aj. This transition effect would 
be ignored if we apply RWA and drop these terms. How- 
ever, as can be seen in the following, the energy difference 
Seij = Ei — Sj is determined by the coupling strengths and 
detunings of the local modes. When the interactions and 
detunings are small enough, 5e,j becomes very small and 
the time dependent factors rotate so slowly that indeed 
these terms cannot be omitted. 
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Notice that /7i„a„ = Ai, U*^Ui„ = UmiU*^ = Smn, we 
can also rewrite the master equation by the local modes 

dtp = t[p,Hs]+CL[p]+CR[p], (6) 

where 

7^ 1 
^Ap] =2-/ ~2^("^^"" ~ -^^o-no-i, p}+) + h.c. 

n 

+ p} + )+h.C., 

n 

and 

i 
i 

Here a takes L,R or 1,3. Cl/r[p] describes the effect 
from the left/right heat bath. For the local modes, be- 
sides the environment induced inter-mode coupling, the 
interaction between two modes can also give the inter- 
mode transition terms straightforwardly in the master 
equation. In this sense, it is concluded that the role of 
the environment mediation is to modify the direct inter- 
mode coupling as a kind of renormalization effect. 

III. LONG-TERM STABILIZATION DYNAMICS 

Comparing with the long time Markovian thermaliza- 
tion process in a heat bath with a single temperature, the 
present complicated environment with two temperatures 
Tl/r cannot stabilize the whole system into a canoni- 
cal thermal state. In this section, we first calculate the 
steady state of the open quantum system by straightfor- 
wardly solving the above master equation Eq. (0]). Then 
we consider the mediating HO as a quantum data bus 
in the large detuning limit. The adiabatic elimination of 
this oscillator can formally induce a direct coupling be- 
tween the left and right HOs. In this case, the analytical 
results about the stabilization can be obtained explicitly. 

A. Steady state in long time limit 

We now consider the indirect coupling case with a me- 
diating data bus. The master equation without RWA can 
be solved with the help of the characteristic function of 
Wigner representation, which is defined as (see Appendix 
B), 

= Tr [p • exp(At -fl-fl^- A)] (7) 
— Tr[p ■ exp(a^ • k — • a)] . 

Here, p, = {ni , p2 , P's)^ and k = {ki, K2, k^)'^ are com- 
plex vectors with respect to the normal and local modes. 



and p = U ■ it. The corresponding Wigner function with 
three modes is defined as the Fourier transform of x(^)j 

With this definition, we obtain the equation of x{p) as 

9tX + z • T • X = z • D • X, (8) 

where z = {pi,p2,P3, pI^P-I^pI), and 



T = 



T ■ 


, D = 





P ' 


T+ 








T and D are 6x6 matrices with 3x3 blocks and P 
defined by 

T^^\{K,-^i)^^^^5^3, (9) 

7;^-^(A7,-A+)+*eA,. 

The equation ([SJ is the Fourier transformation of the 
Fokker-Planck equation about the Wigner function [l^l . 
Formally we give the analytical solution for the steady 
state in Appendix C. Its expression is given as, 

X{p) = xiU ■ k) = exp [z^-i • D' • [zV-^f] , (10) 

where V diagonalizes the matrix T, i.e., V ■ T ■ V^^ = 
diag{di, • • • , dg}, and = [VTtV^],, / {A, + d^). 

All the steady state properties of the composite system 
can be obtained from this formal solution Eq. (|10p . Spe- 
cially, we are interested in the steady state of the HOs 
at the two ends. We can obtain Xa{i^<j) for each local 
oscillator just by setting Hi —Q for al\i ^ a. Notice that 
and D' in the exponent of Eq. ^TU\\ are block diago- 
nal and anti-diagonal respectively, thus it can be verified 
that Xcr{ncr) is always of Gaussian form. Therefore, each 
oscillator finally stays at a state of the canonical distri- 
bution (see the proof in Appendix B), 

x(«:.,<)-exp[-(Arf + (11) 

Since each HO can finally reach a steady state with a 
canonical form, we can treat it as an equivalent thermal 
state and define an effective temperature from its average 
occupation N^^ = {d\da), 

Tf ^uj,/Hl + ^). (12) 

In Fig. 2(a, b) we show the effective temperatures of 
the oscillators at the two ends calculated from Eq. (|10p , 
changing with the detuning IS. — uj^ — ujr and the cou- 
pling strengths ql/r- Here we set ZU = [ujl + '^_r)/2 = 1 
as the energy scale and Wm = 2. 



-0.03 -0.02 -0.01 0.01 0.02 0.03 
Detuning A — ui — lo^ 



-0.03 -0.02 -0.01 0.01 0.02 0.03 
Dctnning A = lji — ujj^ 



-0.03 -0.02 -0.01 0.01 0.02 0.03 
Detuning A — w/, — ujf{ 



FIG. 2: (Color online) The effective temperatures r£ (lower blue ones) and (upper red ones) of the oscillators at the 
two ends calculated from the result without (a) (b) and with (c) RWA. We set cU — {ujl + ijJr)/^ = 1 as the unit, and 
Tl ^ 1,Tr = 3, 7L = 0.002, 7fl = 0.003, LJm = 2. We set ql = Qr = g in (a) (c), and ql ^ 9, Qr ^ 0.8g in (b). We plot four 
groups of curves according to g = 0.02, 0.04, 0.06, 0.08, distributed from outside to inside in these figures. At the small regime 
around the degeneracy point, RWA cannot give us a good enough result, especially when g is weak. The extremum points are 
shifted aside when gL gn- 



When the detuning A becomes large or when their 
coupling strength gL/R becomes small, the two oscilla- 
tors tend to be thermalized with their own heat bath re- 
spectively. The effective temperatures get to the closest 
point around the resonance regime A ~ 0. This observa- 
tion means that they are affected by the reservoir at the 
opposite side and heat transfer happens greatly. When 
gL = gR, the extremum points locate exactly at ujl = ujr, 
while they shift aside when g^ ^ gR- When the inter- 
action becomes strong, the effective temperatures of the 
two oscillators tend to get closer and closer, away from 
that of each heat bath. 

As comparison, we also show some counter-intuitive 
observation resulting from the improper omission of the 
transition terms like 2AipA^- — {A^-Ai, p\ with i ^ j [see 
Fig. 2(c)]. In the large detuning area, this approximation 
shows well consistence with previous result in Fig. 2(a). 
But the effective temperatures always equals at the de- 
generacy point even when the coupling strength g is quite 
weak, i.e., when the oscillators tend to be decoupled from 
each other. Similar counter-intuitive result was also given 
in Ref. Q for a coupled two-spin system. 

With the above comparison, we look back at the 
master equation carefully, the transition terms like 
g-i5eijt^2i,pi] - {A]A,,p}+) with i ^ j contribute to 
the energy exchange of different modes Ai and Aj , which 
can be characterized by (AlAj). The oscillating factor 
at the front describes the phase of this transition. If the 
energy difference Scij = Si — Sj is large, this exchange 
oscillate so fast that this effect vanishes averagely. How- 
ever, indeed Seij is determined by the detuning and cou- 
pling strength. When A and gL,R are quite small, these 
terms cannot be omitted. Therefore, the transition be- 
tween different normal modes plays an important role in 
such a composite non-equilibrium system with Marko- 
vian approximation. The rotating terms take effect to 
the long-term steady state behaviour and cannot be ne- 
glected, especially when the couplings between the sub- 
systems are not quite strong. 

It could be temporally concluded that people should be 



much careful to ignore the transition terms for a multi- 
mode open quantum system in a complicated environ- 
ment. Otherwise, the correct physical effect may also be 
ignored. 



B. Effective coupling in adiabatic limit 

When the detuning of to ujl/r is large, we can ap- 
ply Frohlich-Nakajima transformaliion |15l - [l7| . and elimi- 
nate the mediating degree of freedom to get the following 
simplified Hamiltonian, which describes a system of two 
directly coupled oscillators, 

Hs = lo'l a^ciL + a^a^ -f g{a\aR + CLaJ^), (13) 



where 



dl.R 



L,R 
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i^L.R + 

t^L.R. — t^m 

1 / ghgR ^ ghgR 



(14) 



We can see that the coupling strengths gL and gR con- 
tribute a correction to the renormalized frequencies lo'^ j^. 
For this simplified Hamiltonian, we can write down the 
analytical expression of eigen frequencies Si and the 
transformation U for the normal modes A±. Denoting 
cj^ = aJ — A/2, uj'j^ = lJ+ a/2, we have 



p —a 



2a/3 

f2 -^32 



25 
A' 



(15) 



It follows from Eq. (fTSj) that the energy difference 5eij — 
2Ag depends on the detuning A and coupling strength g. 
When A and g are small, the factors exp[±2iAgt] of the 
transition terms between the two normal modes oscillate 
quite slowly. Therefore, these transition terms may not 
be omitted. 



We carry out the similar calculation for this simplified 
two oscillators system as previously, which gives an ex- 
plicit expression for the steady state of each oscillator, 
described by the characteristic function Xtj^Ka-), 

X.K)=exp[-(7Vf + (16) 
Nf ^[A,NL{e-) + B^NL{e+) 

Here N^^ is the occupation number of the effective ther- 
mal distribution {a — L, R) determined by the linear 
combination of -/Vl.7?(£±), and the coefficients are, 

$ = lLlB.{lL+lBf + l6^l{a^lL+P^lB){l3^lL+a'^'lR), 





1 
T- 



1 
-T 



FIG. 3: (Color online) Two remotely located NAMRs are 
indirectly coupled via a superconducting TLR. The bias volt- 
age provides a difference between the NAMR and the TLR, 
so they can be coupled by a capacitance. At the same time, 
the voltage noise provides each NAMR with independent heat 
bath which has different effective temperature. 



Al = [iLlFiilL + IrT + 16a^Ag(^^7L -I- a^7fl,)7L] , 
Bl = [iLlRilL + IR? + m'^Klia^-iL + P'^Ir^l] , 
Cl = ma^f3^Al{P^-,L+a^lR)lR, 
Dl = Wa^P^Alia^jL + l3^lRhR, 
and 



From the discussion above, we conclude that the inter- 
mode transitions have long-term effect on the compos- 
ite system in a complicated nonequilibrium environment. 
This is quite different from some previous viewpoint 
that the inter-mode transitions only have transient ef- 
fect within the time scale determined by the time-energy 
uncertainty, which usually applies in conventional ther- 
malization process [l|, l9l-[l2| . 



Bn - ma^(3^Al{a^^L+f3^lRhL, 

Cr = [jLlRijL + 7fi)' + l6l3^Al{a^-fL + /3'7i?)7i?] , 

D_R = [jl7r{il + + 16a^Ag(^^7L -I- o^jr^^r] . 



From Eq. (|16p we see that each oscillator achieves a 
canonical state. Especially, at the degeneracy point = 
w'j^, we have = — 1/2, and the difference of the 
populations is. 



eff 



N' 



off 



[NL{e+) - Nn{e+)] + [NLje.) - Nnje^)] 
2(1 + V/7i7fl) 



(17) 



The above equation (jl7p obviously shows that the effec- 
tive temperatures of the two oscillators are not equal at 
the degeneracy point. 

This analytical result Eq. (jl6p gives the same profile 
with Fig. 2 (a, b). From the above arguments, we clearly 
see that the effective temperature of each HO is deter- 
mined by the competition between the interaction of 
the oscillators and that with their reservoirs, which is 
characterized by their ratio /^lIr- When the inter- 
action strength g of the oscillators is large, the differ- 
ence of the temperatures tend to be zero, i.e., the in- 
teraction between the oscillators tend to bind them as 
a whole to be stabilized, and the two oscillators tend to 
get equal temperatures. Moreover, since ^ is renor- 
malized by the interaction energy gL,R, at the degeneracy 
point = LOR, we may have ^ if <?l gR [see 
Eg. (|14l) ]. Therefore, the extremum points are shifted 
aside from ujl = lor as shown in Fig. 2(b). 



IV. REMARKS AND CONCLUSION 

Before we summarize the stabilization result of the 
coupled oscillator system in a non-equilibrium environ- 
ment, we discuss a possible implementation by two nano- 
mechanical oscillators (NAMR) connected via a super- 
conducting transmission line (TLR) to test the theoreti- 
cal results here, as shown in Fig. 3. 

The electromagnetic field inside the TLR may be 
treated as several boson modes whose frequencies are dis- 
cretely distributed [l^,[l^. Only one of the TLR modes, 
which is nearly resonant to the NAMRs, can couple with 
the NAMRs effectively. For example, the voltage distri- 
bution of the lowest even mode along the TLR is. 



V{x) 



27rx c-^. 
cos——{b + ¥), 



(18) 



where — 27r /L\/lc is the frequency of this mode, L is 
the length of the TLR, and Z, c are the inductance and 
capacitance per unit length. The voltage gets maximum 
at the two ends, where the NAMRs are coupled with the 
TLR via a displacement dependent capacitance [2^-[23j. 
The vibration mode of each NAMR may also be treated 
as a single boson. To the lowest order, Cx depends lin- 
early on the movement of the NAMR, Cx — Cxil+x/do). 
Applying a voltage bias Vg to the NAMRs, we have the 
interaction as 

ffint = ^C,"(l + ^)(n^)-K,)'. (19) 

Quantizing the coordinate of the NAMR as a; = Sxo{a + 
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a'''), we obtain an interaction term as Hint = g{a+a^){b+ 
6t). 

The quality factor of the TLR is usually much higher 
than that of the NAMRs, thus we neglect the dissipation 
of the TLR. The voltage noises on the two sides, which 
arise from resistance and bring in the Joule heat, may 
provide the NAMRs with independent heat baths with 
different effective temperatures , and this can be con- 
trolled and measured in experiment [25l - [27| . Typically, 
the frequencies of the NAMRs and TLR are around hun- 
dreds of MHz to several GHz, depending on the sizes of 
the devices, and the coupling strengths are usually sev- 
eral MHz [iMEjIllI- In this regime, it is appropriate to 
make Jaynes-Cummings approximation to obtain linear 
coupling as the above analysis. 

In summary, we have studied the long-term behaviour 
of a coupled HO system connecting with a complicated 
environment which consists of two independent heat 
baths with different temperatures. We deduced a master 
equation with respect to the normal modes. With the 
help of the characteristic description of Wigner function, 
we obtained the numerical and analytical results for the 
steady state of each local oscillator. These results show 
that much attention must be payed when we apply RWA 
in a nonequilibrium system. When the coupling of the 
subsystems is weak, the inter-mode transitions belong 
the RWA may have long-term effect to the steady state 
thus cannot be neglected. This is different from the case 
in conventional thermalization problems, where only one 
canonical heat bath is involved. The non-thermal sta- 
bilization process is determined by the competition be- 
tween the rate of the inter-mode transition and that of 
the energy exchange with each private heat bath. 

This work is supported by National Natural Sci- 
ence Foundation of China under Grants Nos. 11121403, 
10935010 and 11074261, and National 973-program 
Grants No. 2012CB922104. 



Appendix A: Derivation of master equation 

We show the derivation of the master equation Eq. Q 
here. In the interaction picture of Hs + Hb, the interac- 
tion with the environment becomes, 

Vjit)^Vj^{t) + Vfit), (Al) 



where 



j 



(A2) 



Here Ai = Uijdj are the normal modes of the interacting 
oscillators system. 



We take Born-Markovian approximation andput these 
interaction terms into the following equation [ll , 



POO 

dtP = - dTTrB[Vi{t),[Vi{t-T),p(t)<E,pB]] (A3) 

^0 

o 

dT Ttb [V/^ (t) , [y/ (t - r) , p(t) ® ps] ] 
drTTB [vPit), [vP{t-T),p{t)^pB]] . 



Here we assume that the state of each bath is a canonical 
thermal one and does not change with time, i.e., ps = 
Pb ® P% Pb « exp[-77g/fcr^] and = aJk^c^Jk^ 
where is the temperature of the left /right heat bath. 
Thus, terms like Tr b[Vj^ it)Vf [t - T)p{t) (g) ps] always 
vanish, because they only contain the first moment of 
each bath. 

The rightside of Eq. (jA3[) contains two integrals of the 
same form. Each integral gives four terms, one of which 
is calculated bellow as an example. 



d\{t)T,{t) ■ p{t) ®pB-a,{t- T)rt (t - t) 



/o 
where 



dTdi(t)p{t)dAt - r)(rt (i - T)TAt)). 



(A4) 



dl{t)p{t)d,{t 't)^J2 A^i, 



Ji£i-ej)t , isjT 



du! Jc{ui)Nrj{ui) e 



(A5) 



Here, Ja{Lo) = X]k„ l5k„ | <5(w-Wk„) is the coupling spec- 
trum, and N„{£i) = [e-x:p{ei/kTA - 1]"^ is the Planck 
distribution with temperature T„ . The integral Eq. (|A4p 
gives 



Y,U,.U*„A\pA,^ 



,i{ei~Ej)t 



7<t(£j) 



N{ej)+i'P / doj 
Jo 



(A6) 



Here, we denote jai^i) — 27rJ(j(ei), which characterizes 
the coupling strength with each heat bath. The principle 
integral contributes to Lamb shift, and we omit this term 
in this paper. 

The physical meaning of Eq. (|A4|) may be understood 
in the following way. At time t — t, the coupled HO 
system emits energy to the environment, and then ab- 
sorbs back at time t. However, the energy exchange with 
the environment during this process is done by the to- 
tal normal modes Ai but not the local modes a^. Thus, 
when the emission and absorption modes are not the 
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same one, there is an oscillating factor exp[i(ej — ej)t] 
left, deij = Ei — Ej characterizes the splitting amplitude 
resulting from the coupling. By the mediation of the en- 
vironment, the different normal modes Ai of the system 
are coupled together. 

Other terms of Eq. (jA3p can be also obtained as above. 
Each of the two integrals gives the following Lindblad-like 
form with an extra oscillating factor, 



E \UiaUl{lAe^)[Na{Ei) + 1] + 7.fe)[iV.fe) + 



For simplicity, we assume 70- (ei) — 70- does not depend 
too much on lo and can be treated as constant. In sum 
of Eqs. (jA31IA7l) . we get the following master equation in 
Schrodinger's picture, and the oscillating factors do not 
appear. 



dtp^i[p, Y,^^A\A,\ 



A 



+ ^(2A\pA,~{A,A\,p}+), 



(A8) 



where 



^%=^U^^lU*,[NL{Ei)+NL{E,)] 

+ ^-jU,:^U*:,[Nr{e,) + NR{Ej)l (A9) 

^^U,lU*,[NL{Ei) + Nl{e,) + 2] 



UaU*:,[NR{E;) + NR{E,) + 2]. 



The terms with i ^ j describes the transition between 
different normal modes. These terms are often omitted 
by RWA. 



Appendix B: Characteristic function of Wigner 
representation 

The Wigner representation often give us great conve- 
nience to study properties of quantum oscillators. It can 
be defined from a characteristic function i29il , 



Xn,{n.n*) = Tr[e'^^'-^'^p\ 



(Bl) 



The Wigner function is defined as the Fourier transform 

of Xw{n, K*), 



W{a,a*) j d^i^e 



1 ,^ ^ — /ta*+K*a 



X»(At,/t*). (B2) 



For a system that consists of two oscillators, the char- 
acteristic function can be defined as, 



Xi2(ki,K2) = Tri2 [f 



P12 



From this definition, immediately we can find that if 
we had known Xi2(kijK2) for the whole system explic- 
itly, it would be quite easy to get the description of the 
subsystems Xi(2), just by setting K2(i) = in Xi2('«i, ^2), 
without having to calculate the reduced density matrix 
of subsystems Pi (2)- This provides a simple method for 
us to study the state of subsystems. 

Besides, for the thermal state of a oscillator pT = 
exp[— -^a^a], the characteristic function is, 

XT{n,K*)^e^p[-{N+^)\Kf]. (B3) 

Here N = [exp{uj/kT) — 1] ^ is the Planck distribution. 

As seen from the definition, Xui('«, «^*) and W{a,a*) 
can be mapped into each other through Fourier trans- 
formation. It is also well known that there is one-to-one 
correspondence between a physical Wigner function and 
a density matrix. Therefore, there is one and only one 
density matrix p decided by a legal Xw{k, k*)- 

Thus, if we have a characteristic function xi^ji^*) 
which has a Gaussian form like Ea. (jB3l) . with N > 0, 
we can always come into the fact that the corresponding 
density matrix is 



n=0 



(B4) 



where /3®*^ comes from N — [exp(/3'''^w) — l] ^. This 
is a canonical state for the oscillator with an effective 
temperature 1/(3°^. A more rigorous proof lies bellow. 
Proof: From the definition of Xw{k, k*), we have 



(B5) 



62 



e"" "{a\p\a) 



If we have a Gaussian formed characteristic function like 
Ea. (IB3p . we can correspondingly get {a\pT\a) by reversed 
transformation of the equation above. 



{a\pT\a) = I e 

TT 



d n -Ka-+.-a ,g^pr ^^^^^|^|2- 



l 



■ cxp 



(B6) 



iV+ 1 ' L N + V 
On the other hand, we can also expand (a|pT|Q^) as, 

{a\pT\a) ^ E ("1"^-) ("T- \Pt\ n) {n\a) (B7) 



El I \ \ K'^ j -lap 
V m!n! 



8 



Comparing with the expansion of Eq. (jB6p , we can get 
the matrix elements of pT, 



(1 + ^)"" 
(m \pt\ n) = S^n- 



(B8) 



Denote exp[/3'^^aj] = 1 + -^, we can see that pT is a canon- 
ical state. ■ 



Appendix C: Steady solution of Fokker-Planck 
equation 

The standard form of Fokker-Planck equation and its 
characteristic equation are as follows, 



92p 



(CI) 



dt 

/(^, t) is the Fourier transformation of P{y, t), 



The equation of /(^, is a first-order quasi-linear partial 
differential one. It can be solved analytically (s^l, and the 



solution is, 



/(^:t) = <i>tee^'S-..)-exp[--^a,, 



60 



2 ^ ~h A-) 



(C2) 



<&(• • • ) is determined according to the initial condition, 
and $(t oo) = 1. 

In our problem, the equation of the characteristic func- 
tion is 



d 



dtX + z-T ■ — = z • D • z^x 



dz 



(C3) 



The only difference with the standard form is that T is 
not diagonal here. We first diagonalize it and make it a 
standard form. Denoting V -T ■ = diag{di, • • • , d„} 
and z — ■ V , we can transform our equation into the 
standard Fokker-Planck form, 

dtX + ^-A-4^X^^-VT>V^ -^X- (C4) 

Now we could write down the steady solution as 

x(f)=exp[2V-i-D'-(zV-i)^], (C5) 



where D',, = 



d, + d. 
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